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Abstract 

The theory of poles of solutions of P-I is equivalent to the Nevan- 
linna problem of constructing a meromorphic function ramified over 
five points - counting multiplicities - and without critical points. We 
construct such meromorphic functions as limit of rational ones. In the 
case of the tritronquee solution they turn out to be Belyi functions. 



1 Intro 



In this paper we address the following problem in the theory of Painleve 
equations: given a solution y of Painleve I (P-I) 

y"(X)=6y 2 (\)-\, XeC , (1.1) 

and a point w € C, characterize and compute all the inverse images 
of w, i.e. y~ 1 (w) C C. For definiteness we call such a problem the 
Inverse Painleve problem. 

More precisely, we restrict to a rather important instance of the 
Inverse Painleve Problem: We study the poles of solutions of P-I. In 
particular we focus on the tritronquee solution, to be introduced below. 

The main idea behind the present paper is the observation |21] that 
theory of poles of solutions of P-I is equivalent to the Nevanlinna prob- 
lem of constructing a meromorphic function ramified over five points - 
counting multiplicities - and without critical points. 

The first half of the paper is devoted to review this correspondence. 
The second half is dedicated to effectively construct infinitely sheeted 
ramified coverings given the branching locus and the monodromy rep- 
resentation. 

The correspondence between poles and the Nevanlinna problem 
builds on a variety of mathematical constructions which will thor- 
oughly discuss later. For the benefit of the reader we outline the basic 
structure here. 

I Painleve-I is the equation of isomonodromic deformation of a 
system of linear differential equations with an irregular singular 
point 0. This means that given a solution y of P-I and a point a £ 



1 The precise form of the equation is unnecessary for the discussion. It is equation (|3.4|l 
below 
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C, then y[a),y'{a) can be read from the solution of a Riemann- 
Hilbert Problem, provided a is not a pole of y |16| . 
Conversely, suppose a to be a pole of y. Then the Laurent series 
of y looks like 

1 , a(A-a) 2 , (^-«) 3 +5(A _ a) 4 + 0((A _ a)5)i (L2) 



(A - a) 2 10 6 

for some b <= C, with all higher coefficients expressible in terms 
of a, b. 

Remarkably, even though the Riemann-Hilbert problem does not 
have any solution, a careful study |21) shows that linear system 
of equations converge, as A — > a, to the cubic oscillator equation 

ip" = (4z 3 - 2az - 28b) ip . (1.3) 

In a sentence, poles of solutions of P-I are described by the cubic 
oscillator 11. 

II Next step is due to Nevanlinna. Given a cubic oscillator (|1.3p , let 
/ = ipi/ipz the ratio of two linearly independent solutions of it. 
Then / is an infinitely-sheeted branched covering of the sphere, 
ramified over five points (counting multiplicities), and without 
any critical point. 

Nevanlinna remarkably showed that the converse is true |25| : if a 
function f has no critical points and it is ramified over five points, 
counting multiplicities, then f solves, up to affine transformation 
of the z plane, the cubic oscillator il. 3}) for some a, b. 
It follows that solutions of Painleve-I are classified by the branch- 
ing locus of such functions (roughly speaking five points in C 
modulo the action of the Mobius group), while poles of a given 
solution are classified by equivalence classes of monodromy rep- 
resentations. 

Remarkably, the tritronquee solution is classified by the most 
degenerate branching locus, when two pairs of branching points 
happen to coincide and the total number of distinct branching 
points is three; see Theorem [5] below. 

Ill The last ingredient of the theory is the construction of infinitely- 
sheeted coverings of the sphere with given branching locus and 
monodromy representation. After our presentation so far, this 
is equivalent to solving the inverse Painleve problem for poles of 
solutions of P-I. 

Again it was Nevanlinna [22] who showed that they can be con- 
structed as limit of finitely-sheeted coverings, i.e. rational func- 
tions, with the same branching locus and increasing degree, such 
that their monodromy representations approximate the given one. 

The new results contained in the present paper stems from try- 
ing to make effective the Nevanlinna construction of infinitely-sheeted 
branched coverings as limit finitely-sheeted one. Here is a summary. 
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I To illustrate the theory with a neat example, we compute the 
sequences of rational approximants of the ratio of two solutions 
of ip"(z) = (z 2 — 2n — l)tp(z),n e N. These rational functions 
are given two different characterizations that allow them to be 
computed exactly. On one hand, we give the exact location of 
the critical points (see Lemma , on the other hand we give a 
recursive formula (|6.71 16.81) . Latter characterization is still con- 
jectural in the general case (see Conjecture [lj, and proven in a 
particular case (see Theorem [TU)) . 

II We address the Riemann problem of computing a rational func- 
tion given its branching locus and the monodromy representation. 
We give an algebraic treatment of that problem (see Section [7]) . 
We show that the sought function can be described in terms of the 
solution of an overdetermined system of linear equations, whose 
coefficients depend on the unknown location of its critical points. 
After that we show that the Riemann problem can be solved by 
varying the coefficients of the linear systems in order to let it 
have a non-trivial solution, d 

III We study the coverings of the sphere related to the poles of the 
tritronquee solution. We show that the relevant monodromies are 
characterised by two integers n, to. We identify asymptotically 
(i.e. for to, n big enough) the location of the pole with indices 
n, to by making a connection between the monodromy represen- 
tation of cubic oscillators and WKB analysis of the same. 
Eventually, we study the rational approximation of the relevant 
coverings; remarkably they are all Belyi functions, hence they 
are defined over some number field. We compute explicitly some 
of them, in close form, or by solving the Riemann problem nu- 
merically. We warn the reader that our study of these rational 
functions is a preliminary investigation, by no means an exhaus- 
tive one. 

The paper is divided into two parts: Section 2-5 are surveys of 
known results giving a thoroughly exposition of the relation among 
poles of Painleve-I and the coverings of the sphere. The remain- 
ing Sections contain new results. Section 2 introduces the Painleve-I 
equation and the tritronquee solution. Section 3 is devoted to poles 
of solutions of P-I and the cubic oscillators. Section 4 and 5 deals 
with branched coverings of the sphere related to anharmonic oscilla- 
tors. In Section 6 we illustrate the theory with a remarkable example: 
We study the rational approximants of the ratio of two solutions of 
il>"(z) = (z 2 — 2n— l)ip(z), n £ N. In Section 7 we address the Riemann 
problem of computing a rational function given its branching locus and 
the monodromy representation. We give an algebraic solution of that 
problem. Section 8 is devoted to constructing rational approximations 
of those cubic oscillators related to poles of the tritronquee solution. 

2 We expect similar ideas to have already appeared in the literature, even though we 
could not find any paper confirming our expectations. 



3 



Acknowledgements The author is since January 2012 a Scholar 
of the Fundacao da Ciencia e Tecnologia, scholarship number SFRH/ 
BPD/75908/2011. The work started in November 2011 while the au- 
thor was a Research Associate at the University of Sydney, School of 
Mathematics and Statistics, partially supported by the Australian Re- 
search Council Discovery grant nos. DP0985615 and DP110102001. 
We acknowledge fruitful discussions with B. Dubrovin, A. Eremenko, 
N. Joshi, Daniele Masoero and K. McLaughlin. The connection be- 
tween monodromy representation and WKB analysis was prompted 
by questions of B. Dubrovin and A. Eremenko. We also thank A. 
Eremenko for kindly reading and correcting a draft of the paper. 

2 Painleve-I and the Tritronquee Solution 

Painleve I satisfies the Painleve property in a strong sense, as for any 
solution the only singularities in the complex plane are poles. More 
precisely, any solution of P-I is a meromorphic function y(X), A £ C of 
order 5/2 [15]. 

Since all solutions of P-I are meromorphic functions, then for every 
y, w the set is countable and it does not have any accumulation 

point. 

The tritronquee solution is a special solution of P-I. It was found 
by Boutroux in his classical work [S], where he showed that it is the 
unique solution with asymptotics 



The tritronquee solution may also be characterised by other analytic 
properties [TS] or by a Riemann-Hilbert problem [16] . We will use 
latter characterization later on. 

Recently, it has been discovered [8], and partially proven [5], that 
the tritronquee solution describes the leading term of the semiclassical 
expansion of solutions of the focusing Nonlinear Schrodinger Equation 
close to the point of catastrophe. 

In [5], the following conjecture about the location of poles of the 
tritronquee solution, known as Dubrovin Conjecture, is stated. 

Conjecture (Dubrovin). If a € C is a pole of the tritronquee solu- 
tion then | arg a | > 4p ■ 

At the time of writing the conjecture is still open. However, a 
similar but simpler problem about the Ablowitz-Segur solutions of 
Painleve-II equation has been recently solved [3]. 

3 P-I and the Cubic Oscillator 

In this section we recall the relation between poles of solution of Painleve 
I and the cubic oscillator 





(3.1) 
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Here we follow closely our PhD thesis . 
We let 

S k 





2nk 




arg A — 











g|,fc€Z B . (3.2) 

We call the k-th Stokes sector. Here, and for the rest of paper, 
Z5 is the group of the integers modulo five. We will often choose as 
representatives of Z5 the numbers —2, — 1, 0, 1, 2. 

For any Stokes sector, there is a unique (up to a multiplicative 
constant) solution of the cubic oscillator that decays exponentially in- 
side Sk- We call such solution the k-th subdominant solution and let 
ipk(z', a, b) denote it. 

The asymptotic behaviour of ipk is known explicitly in a bigger 
sector of the complex plane, namely Sk-i U Sk U Sk+i- 

Urn ,--H,V e >0. 

2->OC _3 f 4 6 , „ 1 1 ' 

| argz _2|i|<^L_ e Z 4 eXp ^-- Z 2 + - Z 2 j 

Here the branch of z 3 is chosen such that ^ is exponentially small in 

Since ^fe-i grows exponentially in S^, then ipk-i and are linearly 
independent. Then {ipk-i, V'fc} is a basis of solutions, whose asymptotic 
behaviours is known in Sk-i U S^. 

Fixed fe* £ Z5, we know the asymptotic behaviour of {ipk*—i, 4>k» } 
only in Sk*~i U Sfc* . If we want to know the asymptotic behaviours of 
this basis in all the complex plane, it is sufficient to know the linear 
transformation from basis {i/}k—i,ipk} to basis {ipk, V'fc+i} f° r an y ^ € 
Z 5 . 

From the asymptotic behaviours, it follows that these changes of 
basis are triangular matrices: for any fc, Vfc— 1 = V'fc+i + CfcV'fc f° r 
some complex number er^, called Stokes multiplier. The quintuplet of 
Stokes multipliers Uk, k G Z5 is called the monodromy data of the cubic 
oscillator. 

It is well-known (see [21]) that the Stokes multipliers satisfy the 
following system of quadratic relations 

- i(T fe+3 (a, b) = I + (7 k (a, b)a k +i(a, b) , Vfc G Z 5 , Va, b G C . (3.3) 

Hence, it turns out that the monodromy data of any cubic oscillator 
is a point of a two-dimensional smooth algebraic subvariety of C 5 cut 
by 13.31 called space of monodromy data, which we denote by M5. 

Poles of Solutions of P-I It is well-known that P-I is the equa- 
tion of isomonodromic deformation of a system of linear differential 
equations with an irregular singular point. The choice of the linear 
equation is not unique, see for example |17) . 

Here we follow |17| and choose the following auxiliary equation 

-^y- = Q(z ]y ,y',\)iP(z),z,y,y',\€C (3.4) 

Q(z;y,y',X) = 4z 3 - 2Az + 2Xy - Ay 3 + y' 2 



z — y 4(z — y)' 
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We call such equation the perturbed cubic oscillator. 

It turns out that one can define subdominant solutions ?/>fc, and 
Stokes multipliers o~k,k G Z5 also for the perturbed cubic oscillator. 
Moreover, also the Stokes multipliers of the perturbed oscillator satisfy 
the system of quadratic relations p.3[) ; hence, the quintuplet of Stokes 
multipliers of any perturbed cubic oscillator is a point of the space of 
monodromy data M5. 

Since P-I is the equation of isomonodromy deformation of the per- 
turbed cubic oscillator H we can define a map M. from the set of so- 
lutions of P-I to the space of monodromy data; fixed a solution y*, 
Ai{y*) is the monodromy data of the perturbed cubic oscillator with 
potential Q(z;y*(X),y* '(A), A), for any regular value A of y* . The 
map M. is a special case of Riemann-Hilbert correspondence |16j . 

In this paper we are mainly interested in studying poles of solutions 
of P-I. To this aim we cannot use directly the perturbed oscillator in 
this study because the potential Q(z;y,y' , A) is not defined at poles, 
i.e. when y = y' = 00. 

However, this problem can be overcome by studying the auxiliary 
equation in the proximity of a pole. In fact the auxiliary problem has 
a well defined limit. 

Lemma 1. Let a be a pole of a fixed solution y*(X) of P-I and let 
tpkiz] X) denote the k-th subdominant solution of the perturbed cubic 
oscillator \3.J$ with potential Q(z;y*(X),y'*(X), A). In the limit z — > a. 
ipk(z; X) converges (uniformly on compacts) to the k-th subdominant 
solution ipk{z; a, b) of the cubic oscillator 

V>" = (4z 3 - 2az - 286) V> • 

Here the parameter b is the coefficient of the (z—a) 4 term in the Laurent 
expansion ofy*: y* = ^^ + ^2?+ ^l +b {X-af + 0{{X-af). 

Proof. See [20] Chapter 4 or [2T]. □ 

We remark that a similar Lemma for poles of solutions of Painleve- 
II was proven in 

As corollaries of the previous Lemma we were able to prove the 
following Theorems, which define precisely the relation between P-I 
and the cubic oscillator. 0. 

Theorem 1. Fix a solution y* and call er^, k € Z5 its Stokes multipli- 
ers: M(y*) = {a*_ 2 ,..., a* 2 ). 

3 Let the parameters y — y(\),y' — dy j^ of the potential Q(z; y, y' , A) be functions of 
A; then y(A) solves P-I if and only if the Stokes multipliers of the perturbed oscillator do 
not depend on A 

4 Even though the statement of Theorems[T][2]already appeared in [7] by D. Chudnovsky 
and G. Chudnovsky, in |21) we gave a (perhaps the first) rigorous proof |21| . Theorem [3] 
can be proven also by other means |18] . 



The point a G C is a pole of y* if and only if there exists b G C 
such that o~k,k G Z5 are Stokes multipliers of the cubic oscillator 

i>" = (4z 3 - 2az - 286) ij> . 

The parameter b turns out to be the coefficient of the (A — a) 4 term in 
the Laurent expansion of y* . 

Theorem 2. Poles of integrale tritronquee are in bisection with cubic 
oscillators such that o~2 = 0-2 = 0. In physical terminology, these 
cubic oscillators are said to satisfy two "quantization conditions ". 

Theorem 3. The Riemann-Hilbert correspondence Ai is bijective. In 
other words, M5 is the moduli space of solutions of P-I. 

The proofs of above Theorems can be found in [21] and in |20j . 

4 Branched Coverings of the Sphere 

Here we introduce some notions about branched coverings of the sphere 
in order to be able to classify poles of solutions of Painleve-I. Since we 
deal also with coverings from non compact surfaces, the material is not 
that standard. Our main sources are [TU], |26| . 

Given two (not necessarily compact) Riemann surfaces M and N, 
we say that ip : M — > N is a branched covering if ip is a non-constant 
holomoprhic map and if there is a finite set B C N such that ip : 
M \ p^ 1 (B) — > N \ B is a topological covering. The smallest of all sets 
B is called the branching locus. The elements of the branching locus 
are called branch points. 

Two maps tp, \ from M to N are said equivalent if there is an 
automorphism a of M such that ip = \ ct. 

In this paper the target surface will always be C, while the domain 
will be either C or C. We recall that the automorphism of C are the 
affine maps, z — > az + b, a 7^ 0, while the automorphisms of the sphere 
are the fractional linear (or Mobius) transformation z — > ffq^, ab—cd ^ 
0. 

Let Bf be the branching locus of /. Chosen a point w of the 
punctured sphere C \ B and a numbering of f^ 1 (w) (a bijection from 
f~ 1 (w) — > Z„, n possibly 00), we obtain a representation p of the fun- 
damental group 7r(C \Bf,w) in the group permutations of n elements, 
S n . This is called the monodromy representation and it is defined up 
to numbering of the set f^ 1 (w), i.e. up to conjugations in S n . 

4.1 Rational functions 

Holomorpic functions from the sphere to itself are the rational functions 
and every rational function defines a branched covering of the sphere. 

Given a rational function /, z* is a critical point of / if the dif- 
ferential of / evaluated at z* is zero. In a neighborhood of a critical 
point z* , choosing appropriate local coordinates, / has the normal form 
f(z) — (z — z*) k , k > 2. The multiplicity of the critical point z* is 
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k — 1. A point b* is a critical value if it is the value of / at one of its 
critical points. The branching locus of / is the totality of its critical 
values. 

Given / we can read the number of critical points and their multi- 
plicities from the monodromy representation. Let 7& be a small circle 
around the critical value b, then its image p("fb) hi S n is the product 
of some distinct cycles. Every non trivial cycle correspond to a critical 
point and the multiplicity of that critical point is equal to the length 
of the corresponding cycle minus one. 

The number of critical points (counting multiplicities) are related to 
the degree of the function by a simple relation called Riemann-Hurwitz 
formula: Let z\ , . . . , z m be the critical points of / and let v\ , . . . , Vk 
the respective multiplicities. Then 

deg/ = i + ^|i^ (4.1) 

Notice that the Riemann-Hurwitz formula is a non-trivial relation in 
the cycle structure of the monodromy representation of /. 

Equivalence classes of rational functions are in bijection with all 
monodromy representations satisfying the Riemann-Hurwitz formula. 
This is the celebrated Riemann Existence Theorem, which we take from 

HI 

Theorem 4 (Riemann's Existence Theorem). Fix a finite set BcC 
and a morphism p : 7r(C \ B, x) — > S n , for some finite n such that 

(i) the image of tt(C \ B, x) is a transitive subgroup of S„ 

(ii) for any b £ B, if ^b € tt(C \ B) is (conjugate to) a small circle 

around b then p(jb) is not the identity. 
(Hi) the monodromy representation satisfy the Riemann-Hurwitz for- 
mula {4- ty - 

There exists a rational map f , unique up to equivalence, whose branch- 
ing locus is B and whose monodromy representation is p. 

Hence equivalence classes of rational functions whose branching lo- 
cus is B are in bijection with conjugacy classes of monodromy repre- 
sentations satisfying (i,ii,iii). 

4.2 Meromorphic functions 

Holomorphic functions from the complex plane to the sphere are called 
meromorphic functions. Notice that not every meromorphic function 
is a branched covering of the sphere according to our definition. 

In what follows by a meromorphic function we always mean a tran- 
scendental one. Following closely [3] we introduce the concept of tran- 
scendental singularity. 

For any a £ C and disc D(r,a) of radius r > centered in o, 
let U(r) be a component of f~ 1 (D(r, a)) chosen in such a way that 
U(ri) C U{r2) for any r\ < V2- Then either (a) or (b) holds: 

(a) n r >ot/(r) =z,z£C and f{z) = a. z is either an ordinary point 
or a critical point of /. 
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(b) n r >o^(f) = 0- Then the map r -S- [/(r) is a transcendental 
singularity of / . We say that the singularity lies over a. A 
simple argument [5] shows that a is an asymptotic value of /, 
i.e. a is the limit of the function / evaluated on a continuous 
curve tending to oo. According to the Iversen classification, the 
transcendental singularity is either 

(bi) direct, if there exists r such that f(z) ^ a for any z G U(r). In 
this case / : U(r) — > D(r, a)\a is the universal cover of D{r, a)\a. 
Or 

(bii) indirect, otherwise. 

Due to Hurwitz Theorem [TU] we have that / : f^ 1 (D) — > D is a 
covering map if and only if D contains no critical value and no asymp- 
totic values. Then, according to our definition, / is a branched cover- 
ing if and only if the union of the critical values and of the asymptotic 
values is a finite set. 

In this class of functions, the most important Theorem. 

Theorem 5 (Nevanlinna, Elfving). Let f be a function with p < oo 
transcendental singularities and m < oo critical points, lying over q > 2 
points. Then all transcendental singularities are direct and 

where r is a rational function of degree less or equal than p — 2 + 2m. 

In particular, suppose the function f has no critical points, i.e. 
m = 0. Then r{z) is a polynomial of degree p — 2. 

Proof. For the case m = see [33] ■ For the general case, see [§]■ □ 

The operator {f(z), z} is called the Schwarzian derivative. Formula 
(|4.2p shed some light on the connection between branched coverings of 
the sphere and linear differential equations of the second order. Indeed, 
it is well-known that a function / satisfy {f(z), z} = k(z) if and only 
if / = yi /j/2 where yi , yi are two linearly independent solutions of 
y"(z) = -2k(z)y(z). 

Conversely if a function satisfy (|4.2|) with k(z) being a rational 
function then / has a finite number of critical points and transcendental 
singularities (all being direct). We show this below in the case, the only 
relevant for us, where k(z) is a polynomial. 

4.3 Anharmonic oscillators 

By anharmonic oscillator we mean a Schrodinger equation with a poly- 
nomial potential V rn of degree m > 1 

y"(z) = V m (z)y(z) . (4.3) 

Equivalently we mean a meromorphic function whose Schwarzian deriva- 
tive is a polynomial. In fact a function satisfy 
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if and only if it is the ratio of two independent solutions of (|4.3p . 
Particular cases of anharmonic oscillators are the Airy equation m = 1, 
the harmonic oscillator m = 2 and the cubic oscillator m = 3. 

Here, and for the rest of the paper if not otherwise stated, we always 
suppose that the polynomial V m is monic and normalized as follows 

V m (z) = z m + a 2 z m - 2 + ■■■ + a rn , m > 1 . (4.5) 

Following [TT], we outline the coverings structure of the general 
solution of the anharmonic oscillator. With this normalization, the 
Stokes sectors are defined as 



S k = A 



2-Kk 

arg A 



<||,fceZ m+2 . (4.6) 

Let / be a solution of (|4.4[) . Then / satisfies several remarkable prop- 
erties 

• / has no critical points. 



Along any ray r = pe l{p , p G M + contained in S 



lim /(r) = w, . (4.7) 



Moreover, {u>o, . . . , u> m +i} are all the asymptotic values of /. 
Hence / has m + 2 transcendental singularities (all direct) and 
no critical points. 

Wj ^ Wj+i for any j G Z m+2 

Due to Theorem [S] and the previous discussion, we have the follow- 
ing Theorem |25| . 

Theorem 6 (Nevanlinna). A meromorphic function g hasm + 2 tran- 
scendental singularities and no critical points if and only g = f(az + b), 
{f(z), z} = V rn (z), for some a,6eC,o^0 and some V m of degree m. 

Thanks to Nevanlinna theory, we are now able to characterize 
poles of the integrale tritronquee in terms of branched coverings of the 
sphere. Indeed from Theorem[5]we know that poles of integral tritron- 
quee are in bijection with cubic oscillators such that 02 = c-2 = 0. 
The relation 02 = cr_2 = is readily translated into a relation among 
the asymptotic values q in fact it is equivalent to require that two 
pairs of asymptotic values coincide, namely w\ = w_2,u>_i = u>2\ here 
the Wj are defined by (|4.7p . Hence, poles of the integrale tritronquee 
are classified by coverings of the sphere with five transcendental singu- 
larities lying over three points and no critical points. More precisely 
we have the following 

Theorem 7. Poles of integrale tritronquee are in one-to-one corre- 
spondence with ( equivalence classes ) of branched coverings of the sphere 
f : C — > C, with five transcendental singularities lying over 0, 1, 00 and 
no critical points. 



5 The general transformation between Stokes multipliers and asymptotic values was 
given in |23| . We will not need it here. 
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If we normalize the covering map f in such a way that {/, z} = 
4z 3 — 2az — 286 and wq = 0, Wi = iu_2 = 1, twa = f-i = oo. Then a is 
a pole of the integrate tritronquee. 

Proof. Theorem [2] and Theorem [6] □ 

In Section[5]below, we will classify the branched coverings related to 
the poles of the tritronquee by exhibiting explicitly their monodromy 
representations. To this aim, we introduce some elements of the com- 
binatorial theory of coverings. 

5 Combinatorics of Coverings 

In this section we introduce some elements of the combinatorial the- 
ory of branched coverings of the sphere. Here we follow mainly |26) . 
[9|. As an application of this theory, we will show (following |25j ) 
how infinitely-sheeted coverings can be constructed as limit of finitely- 
sheeted ones. 

In this section, if not otherwise stated, we suppose that the branch- 
ing locus bi , . . . , b n is given and ordered and that we have fixed an 
oriented Jordan curve 7 passing through all the branching points, re- 
specting the ordering. This curve divides the sphere into two polygons, 
the inner and outer polygon with sides the arcs (61, 62), ■ • ■ (b n , b\). 

Let / have branching points b\, . . . , b n) we call Critical Graph of 
/ the inverse image of the curve 7 and of the branch points bi , . . . b n 

The Critical Graph allows us to read the critical structure of /. 
In fact, a critical point of / is a vertex of the critical graph whose 
valence is bigger than 2. The multiplicity of a critical point is half the 
number of those faces minus one. The valence of the vertex is also 
the number of the distinct polygons meeting at the vertex. Similarly, 
a direct transcendental singularities is an asymptotic direction of the 
complex plane where infinitely many polygons meet: If we map the 
plane to the unit disc, they become vertex lying on the unit circle 
whose valence is infinite. 

In those cases where the critical points lie on the circle of radius 1, 
then we always choose 7 to be the circle of radius one: The Critical 
Graph is simply the union of the level curves |/(z)| 2 = 1. 

The dual of such a picture is called line complex |26) . Choose a 
point in the inner polygon Pi and a point P in the outer polygon, and 
for any side (bk,bk+x) a line i& connecting Pi and P a crossing just the 
side (pk,k k + 1) and only once (such a line is unique up to homotopy 
in C\{6i, . . . , b n }). The line complex is the inverse image of the union 
of the lines If.. The vertices are bipartite being either inner (inverse 
images of Pi) or outer (inverse images of P ). 

From the line complex we can easily compute the monodromy rep- 
resentation. For any k, we define a map from the inner vertices to the 
outer vertices (and its inverse): //fc(i>) is the vertex adjacent to v w.r.t. 
the edge I}.. The composition Hk{v) = fJ-k f-k-iiv) of two maps is 
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a permutation of the inner (or of the outer) vertices. It defines the 
monodromy representation of a circle around the branch point bk and 
base point Pi (or P a ). 

Notice that the line complex of / (up to orientation-preserving 
homeomorphism of the domain) is far from being unique if n > 2, 
as it depends on the cyclic ordering of the branch points that we have 
chosen. In fact the braid group acts transitively on the possible line 
complexes of a function; see [H] for the case of rational functions, and 
for the case of meromorphic functions. 

An Example We consider the rational function / = 1 3 + 3 3 z 2 ~+f z -i~+z* • 
It has three critical points 1, i, — i with multiplicity 2,1,1 and critical 
value Below we show on the left its line complex and on the 

right its Critical Graph. Since the critical values lie on the circle of 
radius 1, we choose the circle to be the Jordan curve passing through 
the critical values. In purple is \f{z)\ < 1, in grey \f{z)\ > 1. 




To the aim of classifying coverings of the sphere, we formalise the 
notion of a line complex (for more details, see |12| ) and we show for any 
line complex on C (resp. C) there is a branched coverings from C (resp. 
C) to C with that diagram. In other words, we can restate Riemann 
Existence Theorem and Theorem [6] in terms of line complexes. 

An abstract line complex of index n > 2 is a connected multi-graph 
embedded in the oriented C (resp. C) satisfying the following rules: 

• The set of vertices V is a finite or countable set and there is no 
accumulation point of vertices. 

• The vertices are bipartite, V = V,Y[V : There is no edge con- 
necting two vertices belonging to the same partition Vi or V 

• Every edge takes a value from {1, ...,n}. The set is given the 
natural cyclic order 1 < 2 . . . n < 1 . 

• Every vertex is the endpoint of n edges, one for each value. 

• The edges around a vertex belonging to Vi have positive circular 
order. The edges around a vertex belonging to V have negative 
circular order. 

• For any k = 1, . . . , n, we define a map from Vi to V a (and its 
inverse). fJ>k(v) is the vertex adjacent to v w.r.t. the edge k. The 
composition £&(«) = fj,k ° °f two maps is a permutation 
of the inner (and of the outer) vertices. We suppose that for any 
k, Sfc is not the identity. 
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We consider line complexes up to (orientation-preserving) homeomor- 
phisms of C (resp. C). 



5.1 Rational functions 

Rational functions are conveniently represented by line complexes on 
C. In fact 

LEMMA. For every abstract line complex on C of index n, there exists 
a rational function f (unique up to equivalence) with branching locus 
B — {b\ 1 . . . , b n } and that diagram. 

Proof. This is essentially equivalent to the Riemann Existence Theo- 
rem. See [9] for more details. □ 

5.2 Anharmonic Oscillators 

We are interested in studying the line complex of a function whose 
Schwarzian derivative is a polynomial (|4.4|) . Due to Section POl such 
functions are defined by having no critical points and m + 2 direct 
transcendental singularities. 

The branching locus of / is the set of all asymptotic values, namely 
{wq, . . . , w m+ i} = {b±, . . . , b n }. The line complex of / is an infinite 
diagram on C satisfying 

1 The index of the diagram is the cardinality of {wi, . . . , «; m +i} = 



2 For any permutation all non trivial cycles are infinite. 

3 Let Efc be the permutation corresponding to the asymptotic value 
u>i = bk, then its cycle decomposition has as many non-trivial cy- 
cles as Stokes sectors with asymptotic value wi : the total number 
of non-trivial cycles is m + 2. 

The following Theorem due to Nevanlinna establishes the bijection 
between line complexes with properties (1,2,3) and anharmonic oscil- 
lators. 

Theorem 8. Given any abstract line complex of index n on C satisfy- 
ing (1,2,3) as above whose total number of non-trivial cycle is m + 2. 
there exists a meromorphic function f , unique up to equivalence, with 
m + 2 transcendental singularities lying over {pi, . . . b n }, without crit- 
ical points and with the given line complex. 



An example As an example, we classify the line complexes corre- 
sponding to harmonic oscillators y"(z) = (z 2 — E)y(z) with just three 
distinct asymptotic values - they are those line complexes of order 3 
with four non-trivial cycles, all of them infinite. 

Lemma 2. The diagrams T> n ,n £ N represented in Figure]]} are all 

the line complexes of harmonic oscillators y"(z) = (z 2 — E)y(z) with 
just three distinct asymptotic values. 



bn} 



Proof. [2>]. 



□ 



13 



xiV 




Figure 1: The line complex T> n . We omit the labeling of the edges. 

5.3 Rational approximations to Anharmonic oscil- 
lators 

We have seen that branched covering of the sphere with a Unite number 
of transcendental singularities and no critical points -i.e. anharmonic 
oscillators- are in bijection with line complexes satisfying properties 
(1,2,3) of Theorem[Hl 

The problem we address here is how to construct a meromorphic 
function / given its line complex. In light of what we have shown so 
far this is equivalent- if we restrict to cubic oscillators- to computing 
the location of poles of solutions of Painleve-I. 

A solution of this problem was given by Nevanlinna [53]. The idea 
is, given the line complex V of /, to construct a sequence of rational 
function fk whose line complex T>k approximates T>. As was already 
said we suppose that the branching locus is fixed and a curve 7 passing 
through it. Also we choose two points Pj and P Q belonging to the inner 
and outer circles cut by the Jordan curve. 

The construction goes as follows: 

(i) Choose a a vertex vq of V and finite portion Bo of T>, to which 
vq belongs. 

(ii) Similarly let E>k the portion of V whose vertices and edges have 
distance not more than k from Bq. 

(iii) Embed Bk in a disk. 

(iv) Reflect Bk on the circle and switch the bipartition of the vertices 
to obtain a graph Bk with the correct orientation of the edges. 

(v) Finally, on C draw the figure Bk + Bk , by connecting every vertex 
in the frontier of Bk with the corresponding vertex in Bk ■ 

The resulting figure is a line complex T>k on C whose Critical graph 
has a marked polygon, the one to which vq belongs. 

Let fk be the rational function with line complex T>k fixed by re- 
quiring: belong to the marked polygon and i)/fe(0) = Pi, ii)/£(0) = 1, 
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while oo is an arbitrary point belonging to mirror image of the marked 
polygon, for example /fe(oo) = P a . 

Theorem 9. The sequence fk constructed above converges, uniformly 
with respect to the spherical metric on compact subsets of C, to the 
function f with line complex T> 

Proof 0. □ 

As an example we construct here explicitly a sequence 2?j? , k £ 
N, related to the line complex T) n of the special harmonic oscillator 
constructed in Lemma [5] and Figure [1] They are shown in Figure [5] 
below. 



6 The Harmonic oscillator 

To illustrate his theory, Nevanlinna constructed explicitly [24] -with 
ingenuity and owing to the big 1, m +2 symmetry of the problem- the 
sequence of rational approximations of anharmonic oscillators whose 
potential is a monomial 

{f,z} = ~2z m ,m e N . 

For the same purpose, we go a step further. Recall from Section 14.31 
that a meromorphic function satisfying 

{f,z} = -2(z 2 -E) (6.1) 

has four (counting multiplicities) asymptotic values Wq,Wi,W2, W—i, at 
least three of them distinct. Here we concentrate on those / such that 
two asymptotic values coincide, namely 

wq = lim f(z) = W2 = lim f{z) (6-2) 

z— f+oo z — y— oo 

Latter condition is equivalent to the eigenvalue problem for the har- 
monic oscillator: 

y"{z) = (z 2 - E)y(z) et 
lim y(z) = lim y(z) = . 

z— >+oo z— >— oo 

From standard Quantum Mechanics we know that this problem has a 
solution if and only if the 'energies' are quantized according to the rule 
£ = 2ri+l,?ieN. 

We choose wq = w-2 = 1, wi = i, W—i = i and call f n (z) the solution 
of ([6. 2ft ; we also choose the unit circle to be the Jordan curve needed 
to construct the line complex. 

With such a choice, the line complex of f„ is the diagram T> n that 
was constructed in Lemma [5] and Figure Q]|f|. Following the general 
procedure (Theorem for any diagram D n we build a sequence of 
approximating diagrams T> n k, k € N on C. These are shown in Figure 

_m 

6 This follows from Sturm-Liouville theory as the reader can verify. 
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Figure 2: The line complex 

Lemma 3. Let f nk denote the (equivalence class of) rational function 
with line complex T> n k ■ Then f n k is a function of degree 2n + 4k + 2 
whose critical data are as follows 

(i) fn,k have 4 critical points, Zq, Z\, Z2, Z-\ and three distinct critical 
values bo = 1, b\ = i, b-% = —i. 

(ii) There are two critical points zq and Z2 with the same multiplicity 
2k and the same critical value bo- And two critical points Zi, Z—\ 
with multiplicity 2n + 2k + 1 and critical values b\ and 6_x- 

(Hi) If we choose zq = 1, Z\ — i, z_i = —i then Z2 = — 1 and 
fnk(l/z) = l/f nk (z) = f nk (-z) 

Proof. Points (i,ii) are evident from the line complex T> n k- To prove 
(iii), we hx zq = 1, Z\ = i, Z—i = —i and notice that f n k{l/z) = 
= fnk(—z), due to the symmetry of the line complex. There- 
fore Z2 = — 1. □ 

Functions /„ /. have just three distinct critical values and therefore 
are called Belyi functions |19| (more about them in Section [7]). They 
can be constructed recursively using a -partially conjectural- recursion 
formula. First we tackle the case n = 0. 

Theorem 10. Let us normalize fo k by choosing zq = 1,Z\ = i,Z2 = 
—l,Z-\ = —i and let /oft = ■q^, Pok,Qok without common factors. 
Then we can normalize Pok,Qok in such a way that they satisfy the 
following differential relations 

P(f k {z) = V k {z)P ok - l {z),Q' \{z) = V k {z)Q ok - 1 {z), 
V k {z) = (4fc + 2)(z 2 (4fc + l)-l) . (6.3) 

Pofc(0) = -Qofe(0) = l 

Ik 



Pok(0) = Qvk(0) = -lk, 



1=1 



4:1 + 2 

41 



(6.4) 
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And hence 

Q>> k (z) _ Pofe-i 



Q'iM Qok-i 

Moreover, as k — > +00, / ofe (^=) converges uniformly on compact 
subset of C to fo(x) - i.e. the solution of 

{/o, x} = ~2(x 2 - 1) , w = w 2 = 1, wi = i, = —i . (6.5) 

Proof. For sake of convenience we define at = — = ft fc T! ■ Since 

. 9fc /ofe+l 

f k = Qk ~ Pk and Qo fc = 9fc + Pfc relations (|6.3p holds for p, q if and 
only if they hold for P, Q. The only effect of the transformation is to 
map the critical value 60 = 1 — > 0. 

As /ofe (Lemma [3] (hi)), also the function g k has some useful sym- 
metries, namely gk(-z) = gk{\) = -9k(z)- 

Since has degree 4fc + 2 and z = ±1 is a zero of order 2k + 1, 
then pfe = c(z 2 — i) 2fe + 1 for some c <G C*. We choose c = 1. Due to the 
symmetry f k (—z) = —f k {z), q k is an odd polynomial, dcgqk < 4fc+ 1. 

By an explicit computation we verify that (|6.3p holds for p^. We 
prove that the same identities hold for the polynomial q^. To this 
aim we consider the polynomial Rk(z) = q k (z) — Vk(z)qk—i(z), whose 
degree is at most 4k — 1. Due to Lemma [3] (hi), 

Pk (z) T iq k {z) = ±1 + 0({z T i) 2k+2 )~ik > 0. (6.6) 

Due to the relation p" = Vk{z)pk-i then (16.61) implies that z = ±/ is a 
zero of order 2k for R k {z). Hence Rk(z) vanishes identically. 

We prove now (|6.4p . We notice that it is equivalent to Pfc(0) = 
— l,p' fc (0) = 0,<7a,(0) = 0,q' k (0) = 7fc. Since the first three rela- 
tions are already proven, we must verify the latter. Let qk(z) = 
a k z ik+1 + Oiiz^" 1 ). Due to (O]) then ^ = Hence a k = 

a onf=i ^4T~- Remarkably, q' k (0) = a k . In fact, due to the symme- 
try g k (l/z) = —g k (z) and since z 4k+2 pk(l/ z) = —pk(z), we have that 
qk{z) — z Ak+2 qk(l/z) => q' k (0) = a k . By a direct computation we 
verify that go{z) — z 2 ~ 1 => oq = 2. Relation (|6.4p is proven. 

We note that due to (|6.3|) if, under the scaling x = 2\fkz 1 p k and 
q k converge to some entire non-zero functions, then their limits must 

satisfy y"(x) = (x 2 — l)y(x). Trivially the polynomials p k converge to 

2 

—e~ z , which is indeed a solution of the Schrodinger equation. Since 

8fc (0) = 0,^(0) = 7 fc = 4 ^r(fcti) ■ we have that «*(°) '^ ) 
Hence if (7 — fc converges it converges to a non-zero functions. 
We can prove its convergence using (|6.3[) . First we change variable 
and choose Sk(t) = —iqk(it/{2\fk)) so that in (|6.3p Vfe becomes (4k + 

2)(z 2 (4fc + 1) + 1). We let s k (t) = ESo & ?* 2 ' +1 - Usin S §M and 
so = 2t, we prove by (a double) induction that there exist Ci,l € N 
such that 

0< < 6f < Ci, VZ,fc 
Gi+i< ^-(Ci + C-i), V/. 
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The first inequality proves that bf converges to some positive value bf 
while the second shows that s k converges to some function with an 
infinite radius of convergence. Hence g k converges in the scaling limit 
to a solution of (|6.5p . We let the reader verify that wo = u>2 = 0, 
W\ = i, W-i = i. □ 

We leave the generalization of Theorem [TU] for the general f nk as a 
conjecture. 

Conjecture 1. Let us normalize f n _ k by choosing zq = 1,Z\ = i,Z2 = 
—l,Z—i = —i and let f Utk = q 1 ^, P n k,Qnk without common factors. 
Then we can normalize P n k , Q n k in such a way that they satisfy the 
following differential relations 

= V nk {z)P nh -i(z) , Q'U{z) = V nk {z)Q nk . 1 {z), 
V nk {z) = (4fc + 2n + 2)(z 2 (4/c + 2n+l)-2n-l). (6.7) 



Pnk(0) = (~l) n+1 Qnk(0) = l 

Ink 



PM = (-l) n Qo fc (0) = -ink 
k 



(2n + 2)1 



1=1 



M + 2n 



(6.8) 



And hence 



k-1 



ik-1 



Moreover, as k — > +oo, /nfe(^^) converges uniformly on compact 
subset of C to f n (x) - i.e. the solution of 



{fn, x } = -2(x 2 - 2n — 1) , w = W2 = 1, wi = i, W-i 



(6.9) 



Notice that (|6.81 16. 9p can be proven, provided (16.71) holds, using 
essentially the same proof of Theorem [TU] Relations (|6.8[) are rather 
interesting because they solve what, in the large k limit, is the central 
connection problem for the harmonic oscillator. 

REMARK. We stress that Lemma \3(i,ii,iii,iv) allow us to compute 
exactly the functions f„ k , without invoking the differential relations 
fi6.!A 16'. 7)) . How this is possible will be shown in the next Section of the 
paper, see Theorem \ll\ 

Here we list few of the functions f nk (z) and we plot the correspond- 
ing Critical Graph, which is the the inverse image of |/„fc| = 1. In the 
images below, in purple is \f nk (z)\ < 1, in grey is \f nk (z)\ > 1. We let 
the reader associate to any image, the corresponding function. 

l-32-3z 2 +2z 3 +3z"-3z 5 -z 6 



/oi(*0 

fn(z) 
hi{z) 
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-l-3z+3z 2 +2z 3 -3z 4 -3z 5 +z 6 

4-15z-20z 2 +20z 3 +40z 4 -58z 5 -40z 6 +20z 7 +20z 8 -15z 9 : -4z 
-4-15z+20z 2 +20z :5 -40z 4 -58z 5 +40z 6 +20z 7 -20z s -15z a +4z 1(l 

3+16z-36z 2 -48z 3 +50z 4 +48z 5 -36z 6 -16z 7 +3z 8 
3-16z-36z 2 +48z 3 +50z 4 -48z 5 -36z 6 + 16z 7 +3z 8 

2-15z-50z 2 + 100z 3 + 140z 4 -154z 5 -140z 6 + 100z 7 +50z 8 -15z fl -2z 10 
-2-15z+50z 2 + 100z ;i -140z 4 -154z 5 + 140z 6 +100z T -50z 8 -15z a +2z ll) 
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7 The Riemann problem 

In the previous Sections, we have shown that computing the location 
of poles of Painleve I transcendents is equivalent to computing some 
infinitely-sheeted branched coverings of the sphere. These in turn can 
be constructed as limit of a sequence of coverings of the sphere with a 
finite number of sheets. In other words, we have reduced the original 
problem about Painleve transcendents to the construction of a rational 
function once its critical values and its monodromy representation are 
given. We name this latter problem the Riemann problem following 
Plemelj's nomenclature. Plemelj solved the Riemann problem by what 
is now called Riemann- Hilbert technique [27 . We give here another 
treatment, which is algebraic in nature. 

To start with, we introduce a simplification of the Riemann prob- 
lem, that may admit more than one solution, but that can be effectively 
solved by purely algebraic means, as we show below. 

Definition. We call Hurwitz problem the problem of constructing a 
rational function once the following data are given: the critical values, 
the number of corresponding critical points and their multiplicities. 

In general - if the number of critical points is larger than three - a 
Hurwitz problem do not uniquely determine a single equivalence class 
of rational functions but a discrete family of them, because a finite set 
of inequivalent monodromy representations may share the same cycles 
structure 0. 

Let us discuss the solution of the Hurwitz problem. 

7 A Hurwitz problem have no solutions at all, if the data of the problem do not corre- 
spond to the cycles structure of any monodromy representation 
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Definition. Let f be a rational function, we call critical data the col- 
lection of all the critical points of f , their multiplicities and their crit- 
ical values. 

Our method of solution of the Hurwitz problem is based on the 
following observation: If we represent the rational function f(z) = 
as a ratio of two polynomials, then the critical data impose some linear 
relation on P and Q, which can be solved to reconstruct /. In fact 

Lemma 4. Given P, Q polynomials with no common factors, let f = 
^,deg/ = n be their ratio. Then 

z* 7^ oo is a critical point of multiplicity v with critical value b* 7^ 00 
if and only if 

P (l \z*) = b*QW(z*), I = 0,...,v (7.1) 
p(«+D( z *) ^ b*Q (v+1 \z*) 

Here P">,Q"' denote the l-th derivative of P and Q. 

z* 7^ 00 is a critical point of multiplicity v with critical value b* = 00 
if and only if 

Q (l) {z*) = 0, l = 0,...,v (7.2) 

q(«+i)(z*) ^ 

Let p(t) = t n P{\),q{t) = t n Q(\). b* = 00 is a critical point of 
multiplicity v with critical value b* 7^ 00 if and only if 

pW(0) - &V°(0), l = 0,...,v (7.3) 

p (i,+1 )(o) ^ &V" +1) (o) 

b* = 00 is a critical point of multiplicity v with critical value b* = 00 
if and only if 

q (l) {z*) = 0, l = 0,...,v (7.4) 

q( v+1 \z*) ^ 

The reader should notice that equations (|7.1[ 17.21 17.31 17.4p are linear 
in the space of ordered pairs of polynomials C[X] C[X]. 

Motivated by the previous Lemma we define a linear operator on 
the space of pairs of polynomials 

Definition 1. Given a triplet (z, v, b) G C X N x C we define a linear 
operator A n (z, v, b) from the space of pairs of polynomials of degree at 
most n to C +1 : 

A„(z, ! /,| ) ):C[I]„fflC[4^C +1 . 

We call A the critical evaluation map and define its action as fol- 
lows: 

• If z ^ 00 and b ^ 00 then A{P,Q) = ©f =0 PM ~{z) - wQ® (z) . 

• If z ^ 00 and b = then A{P,Q) = ®i =0 (~Q {l) (z)) . 
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• If z — oo and b ^ oo, then A{P,Q) = ® v l= oP {l) (0) - wg (0 (0), 
w/iere p(i) = t n P{\), q(t) = t n Q{\). 

• If z = oo and 6 = oo, £/ien A(P, Q) = ©^ =0 (-g (i) (0)), w/iere 
p(i)=t n P(£),g(t) = t n Q(f). 

More generally, we extend the critical evaluation map to sequences of 
triplets {(Zi,Ui,Wi)}i—i t ___ jS : it is the direct sum of the critical evalua- 
tion maps ®iA n (zi,fi,Wi) 

Due to Lemma01 if a sequence {(zi, fj, 6j)} represents the critical 
data of a rational function / = P/Q then (P, Q) belongs to the kernel 
of the critical evaluation map A({(zt, v>i, 

We show that the converse is true in Theorem [TT1 below. 

Theorem 11. Given a sequence {(zi, fi, 6i)}|_ 1 smc/i that 
(i) ${h,...,b s }>3 

(n,) Ei is even 
(wj /or a// i, Efe J = 6j fj < Ei f 
(iwj a// £/ie Zi 's are distinct 

Letn = 1 + ^| — and .4 = i/j, &i)}f =1 ) oe t/ie critical evaluation 

map of the sequence, acting on the space of pairs of polynomials of 
degree at most n. 
Then 

• (a) If (P,Q) is a non-trivial element of the kernel of A, then 
{P,Q) have no common factors and the sequence {(zi, 6j)}| =1 
is the critical data of f = ^ . 

• (b) The kernel of A has at most one dimension. 

Proof. To avoid considering particular sub-cases of the same proof and 
without loss of generality we may suppose zi ^ oo,Vi. 

a) Let (P, Q) belong to the kernel of A. Uniquely up to a norming 
constant, we can write 

s 

(P(z), Q{z)) = Y[(z - Zi ) k 'r{z)(p(z), q(z j) , 
i=l 

where fc, < Vi, r{zi) ^ 0, Vz and p, g have no common factors. 

Suppose p/q is not a constant function: neither p nor g are zero 
and the they are not both constant polynomials. Then f{zi) = bi and 
Zi has multiplicity /i, > Vi — fc, > 0. 

Hence, due to Riemann-Hurwitz, 

At the same time 

dcg / < n — k — deg r , 
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since deg P < n — J^i &i ~ deg r and dcg Q < n — fcj — deg r. 

Hence, for both inequalities to be satisfied, it is necessary that 
ki = 0,Vi and r = const. Hence (P,Q) have no common factors and 
{(zi,Vi t Wi)~}* =1 represents the critical data of / = ^. 

To conclude the proof we must show that p/q is not a constant. We 
argue by contradiction. 

Suppose p/q is a constant; by a Mobius transformation, we can 
always reduce to the case it is 0: We can suppose that p = 0. Then 
p = P = and Q satisfies (^) = 0, i =, 0, . . . , Wj for all i s.t. 6; 7^ 0. 
Hence Q has at least Ylb jtoi^ + ^) zeroes - 

By hypothesis (i,iii), the number of the zeroes is greater or equal 

than ^i, — +2 = n + 1. Hence Q is zero. Hence (P,Q) is the zero 
vector. Which is a contradiction. 

ii)Let (P, Q), (p, q) two elements of the kernel of A n - To prove the 
thesis we must show that / = P/Q is equal to g = p/q. 

Suppose / - j ^ then Zi is a zero of order + 1 of / — g for 
every i = 1, . . . ,p. Hence / — g has at least ^ vi + s = 2n + s — 2 
zeroes, counting multiplicities. However deg / — g < In which is less 
than the number of zeros of / — g, by the hypothesis on s. This is a 
contradiction. □ 



7.1 The solution of the Hurwitz problem 

Definition. Let {{vi 1 6j)}| =1 be a sequence of multiplicities and branch 
points. We call {(fj, 6j)}f =1 an admissible for a Hurwitz problem pro- 
vided 

(i) ft{6i,...,6 s }>3. 

(ii) Vi is even. 

(Hi) for all i, Eb J=6l v i ^ Ei f ■ 

Given admissible data {(Vi, &i)}| =1 , we can solve the Hurwitz prob- 
lem using the tools developed so far, in particular Lemma H] and The- 
orem [Tl] above. In fact, they show that it is sufficient to find all the 
sequences {{zi, Vi, &i)}? =1 such that the associated critical evaluation 

map -restricted to pairs of polynomials of degree n = 1 + ^3"* - has a 
non trivial kernel. 

To this aim we start by normalizing the problem by fixing three 
of the Zj's, say z\,Z2,Zs to be 0, l,oo and we consider the remaining 
critical points as the unknowns of the problem. 

Since critical evaluation map sends 

C[X] n © C[X] n = C 2 "+ 2 -> C', I = + 1) = 2n + s - 2 , 

i 

we have two cases, s = 3 and s > 3. In the first case, the critical eval- 
uation map has certainly a non-trivial kernel and the unique solution 
of the Hurwitz problem is readily solved by computing a vector of the 
kernel. 
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Conversely, if s > 3, we fix a basis of C[A"]„ © CLY]„ and let 
vi, . . . V2n+2 denote the columns of the critical evaluation map matrix 

- notice that i> 2 ; are polynomial functions of Z4, . . . ,z s . Then the criti- 
cal evaluation map has a non trivial kernel if and only if the following 
system of polynomial equations are satisfied: 

* (ui A v 2 ■ ■ ■ A v 2 „+2) = , (7.5) 

*, A being respectively the Hodge dual and exterior product. 

Definition. We say that a solution of i 7. 5\) is admissible if all the 
Zj,j = 1, . . . , s are pairwise distinct. 

Lemma 5. Let a Hurwitz problem with admissible data {{vi, Wi)}f =1 
be given. Then the critical points of a function satisfying the Hur- 
witz problem are admissible solutions of 17.5)) . Conversely for every 
admissible solution z\,...z s of i 7. 5}) the sequence {zi, Vi, &i}j=i,... lS is 
the critical data of a unique function, and this function satisfies the 
Hurwitz problem. 

Proof. It follows immediately from Lemma 0] and Theorem 1111 □ 

Summing up, we have shown that the Hurwitz problem is essentially 
equivalent to the system of polynomials equation (|7.5[) . 

REMARK. A finite set of different admissible solutions of 17. 5\) may 
lead to the same solution of the Hurwitz problem, if the latter has some 
symmetries. 

REMARK. Note that d 7. 5\) is a system of ( 2r / +2 ) equations, where 
I = Yliiyi + l) = 2n — 2 + k in k — 3 unknowns. However, we know 
that the system has only a finite number of admissible solutions since 
the Hurwitz problem has a finite number of solutions. 

The solution of the Riemann problem Given a Riemann 
problem, i.e. a branching locus and a monodromy representation, we 
associate to it a Hurwitz problem by extracting from the monodromy 
representation its cycles structure. The solution of the Riemann prob- 
lem is one of the admissible solutions of the associated Hurwitz prob- 
lem. To select the correct one it is sufficient to compute the monodromy 

- for example by plotting the Critical Graph - of all solutions of the 
Hurwitz problem. However, it may happen that further insights (e.g. 
inequalities on the location of the critical points) about the specific 
problem allow to a-priori select the correct solution without comput- 
ing the monodromy. This is the case, for example, of some Riemann 
problems associated with poles of the integrale tritronquee, see Section 
11] below. 

Belyi Functions We end our discussion of the Riemann problem 
by introducing the Belyi functions. They are defined to be those ra- 
tional functions (or more generally, those meromorphic functions on 
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a compact Riemann surface) with just three critical values, namely 
0, 1, oo. 

The importance of Belyi function lies in the following important 
property. 

Theorem 12. If a rational function f has critical values 0, l,oo, then 
f is equivalent to a function g defined over a number field K , Q C 
A'cQ (here Q is the algebraic closure ofQ). 

In other words, Va; € K , g(x) € K. Equivalently let g = P/Q, 
then we can choose P = ^2a n z n ,Q = ^2b n z n in such a way that 
Vn, a„, b n e K. 

Proof See 0. □ 

Our solution of the Riemann problem give a hint of why the The- 
orem holds. In fact, system (|7.5[) has integer coefficients so that we 
can expect its solutions to be algebraic numbers. However, a full proof 
would require to identify the dimensional component of solutions of 
(|7.5p as well as deep instruments from Elimination Theory, such as the 
Shape Lemma. We will not go into that. 

8 Poles of the Tritronquee Solution 

According to Theorem [3 poles of the Tritronquee solution are in cor- 
respondence with those cubic polynomials V(X; a, b) = 4A 3 — 2aA— 286 
such that the solution of the Schwarzian equation has just three dis- 
tinct asymptotic values, namely 

{/,*} = -2V(X;a,b), 
Wi = U>_ 2 , W-X = W2 ■ (8.1) 

Our aim is to construct such transcendentals functions as limit of ra- 
tional functions. Since the number of distinct asymptotic values, three 
in this case, is the number of critical values of the rational approxi- 
mates, then all the rational functions that we are looking for are Belyi 
functions. 

Contrary to the case of the harmonic oscillator, we will be unable to 
give explicit general formulas as we expect the approximating functions 
to be defined -generically- over some number field different from Q. 

We proceed by classifying all line complexes corresponding to the 
poles of the tritronquee, hence by classifying poles themselves. 

Lemma. The line complexes of cubic oscillators with only three dis- 
tinct asymptotic values are classified by two integers m, n G N X N. We 
let X> n,m denote them. They are represented in Figure 

8.1 Lines complexes and WKB approximation 

Before tackling the construction of the rational approximants. we an- 
swer here the following question: where is, approximately, the pole of 
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Figure 3: The line complex T> n ' 
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the tritronquee solution corresponding to the line complex T> n,m l To 
this aim we have to analyze the connection between the line complexes 
and the WKB analysis of solutions of the cubic oscillators. Here we 
simply state this connection without entering into the details of the 
proof, as it would require the introduction of heavy machinery from 
WKB analysis of our previous works. In this Subsection we suppose 
the reader to be be familiar with the results of ]21^ , J22j/ . 

In [m we have derived a system of two equations (equation (25) of 
|21|). called Bohr-Sommerfeld-Boutroux system, for the unknowns a, b 
coefficients of the cubic potential, describing the asymptotic location of 
poles of the tritronquee solution. Solutions of that systems are discrete 
and classified by two positive integers I, k. 

Given two quantum numbers I, k and the unique solution a l ' k , f3 l ' k 
of the Bohr-Sommerfeld-Boutroux equations, for I, k big enough there 
exists a unique pole a of the tritronquee close (in some precise sense) 
to a l ' k , see Theorem 1 in [2"2"] , 

The answer to the above question is as follows: the pole of the 
tritronquee solution corresponding to the line complex T> n ' m is, at least 
for large n, m, the pole close to the solution of the Bohr-Sommerfeld- 
Boutroux system with indices I = n + l,k = m + l. 

Let us give some hints of the proof. Let us normalize f n > m by 
choosing wq = 0,w± = 1, W2 = oo: then the Critical Graph is the union 
of level curves Ref n,m = and the numbers n, m count the distinct 
curves Ref n,m = that asymptotically lie in the Stokes sectors Sq and 
5*2, S-2 (Stokes sectors as defined in (|3.2p ). 

The WKB approximation of the cubic oscillator [5T], allow us to 
approximate the solution / of {/, z} = —2V(z;a,b) (fixed by choos- 
ing w = 0, w i = l,w 2 = oo), as f(z) - f(z) = e s( - z \S(z) 
i ^V(u; a, b)du, where \/V is defined in the plane with some cuts. 
Clearly given a representation of / as above, the lines Ref = corre- 
spond to ImS = dit, d £ Z. 

Remarkably, the indices Z, k of the Bohr-Sommerfeld-Boutroux sys- 
tem count the number of lines ImS = dir that connect So with S%, 
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s- 2 E 

Hence we have two different counts for the number of of lines Ref = 
connecting Sq with 52 and S-2- We have the theoretical (exact) 
one given by the line complex, and the approximate one given by the 
WKB analysis. For n,m big enough, the two counts can be proven to 
coincide. 

8.2 Rational Approximants 

In the present Section, we start studying rational approximation to 
18.11 As it as already remarked, this is a preliminary study. 

To construct rational approximants to the function with lines com- 
plex X> n m , we follow the general procedure (see Theorem [5]) and for 
any diagram T> n m we build a sequence of approximating diagrams 
V^ m , fceNonC. These are shown in Figure gj We let /™' ln denote 
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Figure 4: The line complex Jf^ rn 
the (equivalence class of) rational function(s) whose line complex is 

Lemma 6. The function f^' m has five critical points zq, Zi, Z2, Z— 1, Z—2, 
three distinct critical values b$,bi = &-2>&-2 = &— 1 and multiplicities 

v = 2(m + n + k + l),i> 1 = 2k, v 2 = 2(n + k) + 1 , 

v-i = 2k,v- 2 = 2(m + k) + 1 . 

Let us fix /I 1 '™ by choosing bo = zq = 00, b 2 = 22 = 0, fo_2 = ^_2 = F 
Then 

(i) fk' m (z) = f^' m (z), where — stands for complex conjugation. 

(ii) fk ,m (— 2 + 1) = —f k ' n (z) + l. Hence, if we denote by z±\(m,n) 
be the critical points Z±\ of the functions f n,m , we have 

Zi(m,n) = z_i(n, m) . (8-2) 

(Hi) Fix wq = zq = 00, W2 — Z2 = 0, W-2 = z-2 = 1 then 

zi,z-t G R, z-i > 1 et zi < . (8.3) 

8 Notice that the Stokes graph of solution of the Bohr-Sommerfeld-Boutroux system 
(Figure 4 in |21| ) is the bifurcation diagram of the lines ReS = 0. 
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Proof. The multiplicities of the critical points as well as the symmetries 
are read directly from the line complex. 

Here we prove (i) and leave (ii) to the reader. The critical graph 
of fj!' m is the level curve Imf^' m . By the construction of the line 
complex, there is a closed union of edges L of the Critical Graph C 
cutting the critical graph in two two halves which are, topologically, 
one the mirror image of the other. Hence C is equivalent to its mirror 
image. 

Let us consider the function g?' = f^ ,m {z). They have the same 
critical values and the critical graph D of g k ' m is the mirror image 
under complex conjugation C. Hence D and C are equivalent and so 
-by the Riemann Existence Theorem- / and g are equivalent. Since 
three critical points Zq, Z2, z-2 are fixed by complex conjugation then 
9k' m = fk ' ■ Moreover, the real axis is the line L cutting C into two 
equivalent halves. 

We now prove (iii). To show that z_i > 1, we look at the restriction 
of the Critical Graph on the line L, i.e. the real axis. Since in the 
Critical Graph there is no edge connecting z-\ and Z2, or connecting 
Z\ and Z-2 then z_i > z_2 = 1 and zi < Z2 = 0. □ 

8.3 

Here we tackle the explicit construction of the functions /"'". All these 
functions are defined on Q. 

Definition. Let A n (z) be the rational functions with critical data 

{(0, 2n, 0), (1, 2n, 1), (oo, 2n, 0)} , n e N . (8.4) 

We call it the n-Airy function. Its line complex is shown in Figure [3 
below. 

Functions A n coincide - up to a fractional linear transformation- 
to the functions fh (case h = n, q = 3) in [24] §2. They are the 
approximants of the solution of the Airy equation 

{/,*} = -2*. 

The coefficients of A n are given explicitely in [24j . 




Figure 5: The line complex of the function A, 
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Theorem 13. /"•"(» = A n {f°'°{z)), / °'° = -2z 3 + 3z 2 . The critical 
data of /"'" are 

{(0,4n+l > 0)(l,4n+l > l),(c»,6n + 2 > oo) > (-i 2n, 1), (|,2n,0)} . 

(8.5) 

Proof. The function ^4 n (/ ' (z)) has critical data (|8.5p so it satisfies 
the same Hurwitz problem as f™' n (see Lemma [B]). To prove the The- 
orem it suffices to show that the two functions have the same mon- 
odromy representation. This can be easily verified by composing the 
line complex of A n and of □ 

The possibility of the non-trivial decomposition /™'"(z) = A„(/q'°(z)) 
is due to the fact that the monodromy group of /™ ,n is imprimitive 
(Ritt's Theorem, see [H]). Whether other functions fu' m have im- 
primitive monodromy group is an interesting combinatorial problem 
that we do not address here. 

Below are the Critical Graph of /"'", n = 0, . . . , 3. Here zq = wq = 
1, z±2 = w±2 = e l= 5~ and the Graph is the level curve |/| = 10. 
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8.4 Numerical Experiments 

We have computed numerically the functions fu' m for small n, m, k, 
by solving the corresponding Riemann problem. To this aim -fixed 
bo = z a = 00,62 = Z2 = 0,fc_2 = 2-2 = 1 - we solved the associated 
Hurwitz problem by looking for admissible solutions of the system 
(|7.5p . For small n,ni,k conditions (|8.3[) are enough to select a unique 

9 In the last picture, 1 is a critical point of order 20, hence close to 1 the function is a 
constant up to machine precision. This is the reason why the plot is not well-resolved. 
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solution of the Hurwitz problem, which is the solution of our Riemann 
problem. 

Once the functions are calculated we normalize them according to 
general prescription (see Theorem [5]) to assure the convergence of the 
k — > oo limit: The sequence f k ' m would converge to the function 
sought function /™< m up to an affine transformation of the domain. 

Due to Theorem [3 for extracting the value of the corresponding 
pole of the tritronquee solution, the function / has to be normalized 
in such a way that 03 = 4, a 2 = and that wo = 1. 

In practice we have computed the zero z* of the second derivative 
of {fk' m , z} closest to z = and applied a dilatation to get {fj!' m , z} = 
-2(4(z - z*) 3 - 2al m (z - z*) - 28&£' m . As k -> 00, a n k ' m -+ a k > m , the 
pole of the integrale tritronquee with line complex 2??' m . 

Numerical computation of /?'° To show how the method works, 
here we present /, ' , k < 5 and we postpone a more systematic pre- 
sentation to a subsequent publication. 

To compute them, we have solved the corresponding Hurwitz prob- 
lem. Due to (|8.2p . system (|7.5p reduces to a system of polynomial 
equations in one variable, the location of the critical point Z—\. 

Below we list the resultant of that system, which is the polynomial 
equation satisfied by Z-\. We call P^ corresponding polynomial. 

P {z) = -3 + 2z, P 1 {z) = -7 + 4z, 
P 2 {z) = 693 - 2310z + 2772z 2 - 1416z 3 + 256z 4 , 
P 3 (z) = 1859- 6006z + 6864z 2 - 3328z 3 + 576z 4 , 
P±(z) = -23056709 + 187166226z - 667288284z 2 + 1371124664z 3 
-1788975552z 4 + 1536508416z 5 - 868265344z 6 
+311072256z 7 - 64052736z 8 + 5767168z 9 , 
P 5 (z) = - 16232365 + 129858920z - 454506220z 2 + 913469840z 3 
-1161895600z 4 + 970042112z 5 - 531587776z 6 
+184365056z 7 - 36708352z 8 + 3194880z 9 . 



All Pk listed above are irreducible over the rational numbers. For both 
P4 and P5 the Galois group is the full Sg. 

Computed z_i, the function f k ' is found by calculating the kernel 
of the critical evaluation map, whose terms are polynomials in z_i with 
integer coefficients. Hence the number field defined by P^ is also the 
field of definition of A' . The approximations a^°,3 < k < 6 of the 
pole a ' are 

0.0 r, r-7 0.0 r, r-o 0.0 n r n 

a 3 =—2,57, a 4 ' = —2,53, a 5 =—2,50. 

The reader will convince herself that a 0,0 is the first pole on the negative 
axis as computed in [TS], namely a ' = —2, 38 . . . I 10 l . 



10 The corresponding solution of the Bohr-Sommerfeld-Boutroux system is a 1 ' 1 — —2.32 
[21 j. This confirms our expectations (see Subsection l8.ll) . 
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We plot the Critical Graph of f k ' ,k — 0, ...,5. Here wq = 
l,w±2 = e l ^~ and the Graph is the level curve |/t'°(^)| = 1. 
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/> ' To give an example of poles that are not real, i.e. n ^ m, we 
plot / fc 1 ' °,fc = 0,l- 
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